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ABSTRACT 



The principal underwater sound energy radiated by terminal velocity raindrops at sea is 
due to micro-bubble entrainment and oscillations which occur for drops of the two di- 
ameter ranges 0.8 to 1.1 mm (Type I) and 2.2 to 4.6 mm (Type II). In the absence of 
bubbles, particularly between 1.1 and 2.2 mm, the impact sound radiation is significant. 
The Type I bubbles radiate at frequencies close to 15 kHz, whereas Type II bubbles ra- 
diate between 2 and 10 kHz, depending on the drop diameter. Therefore Type II bub- 
bles, w r hich are common in moderate to heavy rainfall, offer the opportunity to 
determine rainfall drop distribution and total rainfall rate by remote underwater listen- 
ing. Type II bubbles radiate more energy when the drop and surface temperatures differ, 
e.g., almost twice as much energy when the drop and surface temperatures dilTer by 10 
0 C. Type II bubbles radiate less energy in saline water, e.g., 45 % as much energy at 
a salinity of 35 ppt as for fresh water. The distinctive sound spectral shape for a par- 
ticular diameter raindrop does not change appreciably with extreme differences of tem- 
perature (0 to 22° C ) or salinity (0 to 35 ppt). It is possible, therefore, to condense the 
data acquired from hundreds of drops in our laboratory' into a single relation which gives 
the average energy radiated by a Type II raindrop as a function of drop volume, tem- 
perature and salinity. Using this relation, we find good agreement between measure- 
ments at sea and the predicted sound spectrum for an assumed reasonable drop size 
distribution. Also, the total rainfall rate and drop size distribution has been calculated 
from sound spectra measured at sea (the inverse problem.) These early successes lay the 
groundwork for real time measurements of total rainfall rate and drop size distributions 
in moderate to heavy rainfalls inferred by remote underwater listening. 
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I. INTRODUCTION 



The oceans comprise over 70 % of the Earth's surface. Yet, although meteorologists 
know a great deal about the characteristics of rainfall on land, very little is known about 
the characteristics of rainfall at sea. Many techniques have been used in recent years, 
with varying degrees of success, in an attempt to estimate both the rainfall rate (and 
drop size distribution) over the open seas. Direct measurement by rain gages on 
floatable platforms are expensive and particularly difficult in heavy sea conditions. In- 
direct measurements such as radar and satellite imagery lack independent verification 
of drop size distribution. Clearly, other methods must be developed to overcome these 
obstacles. 

Since Knudsen first measured underwater ambient sound in 1948, precipitation has 
been known to contribute significantly to the overall sea noise spectra. Wenz (1962) 
confirmed these measurements and presented a discussion of the possible noise sources. 
His results, shown in Figure 1 on page 2, show that wind and precipitation noise are the 
primary contributors of noise from 1 kHz to 50 kHz. It had been hoped that precipi- 
tation noise could be separated from wind noise, as their respective spectra have different 
shapes (Lemon and Farmer, 1984), and that a direct correlation between sea noise and 
rainfall rate would occur for this frequency band. However, early attempts to secure 
these results have been only marginally successful (Scrimger,et al., 1989; Tan, 1990). 
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Figure 1. Ambient Noise Sources in the Ocean (Wenz, 1962) 



Franz ( 1959) identified two sources of sound from an individual raindrop: the impact 
of the drop on the water surface, and a bubble, which is sometimes formed after a delay 
of at least 30 msec. This sequence is shown in Figure 2 for a 4.2 mm diameter raindrop. 
The bubble (when it occurs) radiates strongly with a frequency given by (Minnaert, 
1933): 



fo = 



1 

2na 0 




( 1 ) 



where p„ and are the local pressure and density, a„ is the equilibrium bubble radius 
and y is the ratio of specific heats ( c p and c„) for air, and is generally assumed to be 1.4. 
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Figure 2. Typical Impact Followed by a Time-delayed Bubble Drop Diameter 4.2 
mm; Sampling Frequency 250 kHz 



Recent studies (Kurgan, 1989; Pumphrey, et al., 1989; Medwin, et al., 1990) have 
been centered around small raindrops ranging in size from 0.8 to 1.1 mm diameter. This 
diameter of drop will produce a bubble 100% of the time when falling from terminal 
velocity and impacting at normal incidence on a smooth surface. Since a bubble is 
known to radiate several orders of magnitude more energy than the impact, this might 
have been an important region for correlating the rainfall rate with the sound produced. 
Unfortunately, several separate phenomena tend to deteriorate the correlation. 

Wind, which is almost always present at sea, imparts a horizontal velocity to the 
small raindrops. This tends to make the drops impact the surface at an oblique angle. 
Studies have shown that the percentage of bubble creation from small (0.8 - 1.1 mm) 
raindrops decreases with increasing angle of incidence: decreasing from 100% at normal 
incidence to 0% at about 25° (Kurgan, 1989). Also, there are other sources of bubble 
noise in the same frequency range of the small bubbles, such as bubbles produced by 
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breaking waves. In short, these effects make it difficult to predict rainfall rate from the 
sound produced by small drops. 

Total rainfall rate (TRR) in mm hr is determined by the drop size distribution, the 
terminal velocity and the diameter of the drops falling (Nystuen and Farmer, 1989): 

TRR = f Z) 3 n(D) Vj(D) dD (2) 

J 0 o 

where n(D) is the drop size distribution [Drops/m J ,’mm diameter increment ], V-,{D) is the 
terminal velocity [m/s] and D is the drop diameter [mm]. Since the equation for terminal 
velocity can be approximated by V T x 4.6 Di mis over the range of 2.2 mm - 5 mm di- 
ameter (Snyder, 1990) and most raindrops in natural rainfall are less than 5 mm, we 
obtain the rainfall rate (RR) in mm hr for drops with diameters between 2.2 - 5 mm us- 
ing the mixed units favored by meteorologists: 

RR„ = 2.8 x 10' 3 7i f 5 n(D) Z)i dD (3) 

*' 2.2 

Total rainfall rate is strongly dependant on the larger drop diameters. The effort to 
correlate rainfall rate with sound levels therefore shifted to larger raindrops beginning 
with Snyder's work. Some typical drop size distributions are shown in Figure 3. They 
were obtained during a convective storm at Clinton Lake, IL, in October, 1982 (courtesy 
Nystuen) using an Distromet, LTD. distrometer with a drop diameter resolution of 
±0.1 mm for diameters between 0.3 - 0.8 mm, ±0.2 mm between 0.8 - l.S mm, ±0.3 mm 
between 1.8 - 3.3 mm, ±0.4 mm between 3.3 - 4.5 mm and ±0.5 mm between 4.5 - 5.6 
mm. 
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0.6 mm/hr 
12 mm/hr 



52 mm/hr 
92 mm/hr 



Figure 3. Typical Drop Size Distributions for Light to Heavy Rainfall (Courtesy 
Nystuen) 



Initial investigation in the NPS laboratory of the sound produced by large, terminal 
velocity raindrops produced startling results (Snyder, 1990). Snyder was able to directly 
correlate bubble frequency with drop diameter for drops between 2.7 and 4.6 mm (5 mm 
is the largest drop normally found in heavy rainfall). With this knowledge, the goal of 
correlating rainfall rate with sound level appears to be feasible. 

The goal of this thesis is to extend the results of Snyder's work to include both the 
energy radiated by large raindrops and the dependence of the energy radiated on the 
ocean parameters of temperature and salinity. With these findings, it will be possible to 
obtain both total rainfall rate and drop size distribution from the underwater sound ra- 
diated by rainfall. In addition, it will also be possible to construct an underwater rainfall 
sound spectrum level given only the drop size distribution. No other method is currently 
able to accomplish either of these tasks successfully. 
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II. EXPERIMENTAL SETUP 



A. ACOUSTICAL EXPERIMENT 
1. Drop Toner 

The Naval Postgraduate School has a facility unique in raindrop sound research: 
a 3 m x 3 m x 26 m vertical utilities shaft with a 1.5 m deep x 1.5 m diameter redwood 
lined anechoic tank at its bottom (Figure 4). The top of th- utilities shaft is fully ac- 
cessible to allow simulated raindrops to fall through the entire 26 m travel. This distance 
is adequate to produce terminal speed in all sizes of drops found in normal rainfall ( < 
4.8 mm diameter) as shown in previous work (Snyder, 1990). This "drop tower" was 
used in all portions of this research, except for the crossover energy experiment for 
intermediate raindrops with diameters between 1.1 to 2.2 mm. 




Figure 4. Diagram of Drop Toner and Redwood Tank 
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2. Anechoic Redwood Tank 

The tank used to conduct the acoustical portions of this work is a cylindrical, 
Winemaker's barrel with dimensions of 1.5 m height and 1.5 m diameter. Although the 
redwood used in the construction is known to be a good absorber of acoustic energy, 
an additional 15 cm of rough redwood was inserted to further mitigate the effects of 
reverberation. This lining was not present for the data used in Snyder's (1990) work, 
but the appropriate correction factors (to account for the additional reverberation) have 
been used for comparison. A minimum of 7 dB reduction in reverberation level at 3 kHz 
and up to 22 dB reduction at 20 kHz was noted when the lining was used. 

The tank was filled with tap water for all portions of this work, with the excep- 
tion of the salinity experiment. To simulate ocean water containing 35 ppt salinity, the 
tank was filled with Sea Salt, manufactured by Lake Products, Co. The temperature of 
the water can be maintained at elevated temperatures (up to 30 °C) by continuously- 
running an attached filter and covering the tank when not in use. 

3. Low Velocity Drop Apparatus 

To obtain the lower (non-terminal) velocities required for the crossover energy 
experiment, a 3 m ladder was used to achieve a variable height above a 20 gallon fish 
tank. The crossover energy experiment did not require any quantitative acoustical 
measurements, so the reverberant nature of a small aquarium did not affect the results. 
The tank also provided an excellent setting for filming the impact sequence of the drop 
at the air/water interface during the high speed filming, discussed later this chapter. 

4. Drop Producing Apparatus 

The early work, based on the same data contained in Snyder's (1990) exper- 
iments, used several means of drop production. An Eppendorf digital pipette (Model 
4710 0.5 - 10.0 /d) was used for the 2.2 - 2.7 mm drop diameter range, with a published 
volume accuracy of ±1%. For drops in the range of 2.7 - 3.6 mm diameter, an Eppendorf 
digital pipette (Model 4710 10 - 100 n\) was used, with a published volume accuracy also 
of ±1%. For drops with a diameter greater than 3.6 mm, individually calibrated glass 
eye droppers were used. The accuracy of the eye droppers was calculated by measuring 
the volume of 50 drops at least 5 times (collected in a graduated cylinder) and was 
measured to be ±5% by volume. The acuracy of the graduated cylinder was ± 0.1 ml. 
For example, a 4.6 mm drop has a volume of 50 jul, therefore 50 drops should have a 
volume of 25 ml. In all cases of volume measurement, the volume of the 50 drops was 
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within 1.2 nil, or 5%. The accuracy was also verified on an Ohaus precision balance to 
be ±5% by mass for each individual drop. 

In the later work, a smaller range of drop sizes led to the automation of the drop 
producing process. A standard medical intra-venous drip bag was used to feed a cali- 
brated glass eye dropper that produced a stream of separated drops with an adjustable 
drop rate. This freed one person from the task of creating hundreds of drops at a sitting 
and minimized the drop strike radius by maintaining a constant "aim." The repeatability 
of the droppers was calculated by measuring the volume of 100 drops at least 5 times 
and was +5% of the expected volume for all drop diameters used. Again, the accuracy 
was verified using an Ohaus precision balance to be ±5% by mass for each individual 
drop. 

5. Hydrophone 

The construction of the hydrophone was similar to that of an LC- 5 and con- 
sisted of two 1/8 inch coaxial cylindrical barium titanate elements. The hydrophone was 
calibrated by both the spherical reciprocity and comparison calibration methods. Its 
response was quite flat (+ 1.5 dB) up to 50 kHz and ±5 dB from 5 - 300 kHz (Figure 5 
on page 9). It was positioned at 15 cm depth for the early work and 6 cm depth for the 
later work. Although the shallower depth included the near-field acoustic and hydraulic 
effects experienced during an impact, the increased signal to noise ratio provided a much 
cleaner signal. Also, the near field effects can be accounted for, as described in Chapter 
IV. 
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Frequency (kHz) 

Figure 5. Hydrophone Calibration Curve 



6. Amplification and Filters - Early Work 

The signal from the hydrophone was amplified by a PAR 113 Low Noise Pre- 
Amplifier with a gain of 2000 for drop diameters greater than 3.4 mm and 10,000 for 
drop diameters less than 3.4 mm. It was then passed through two Krohn-Hite band pass 
filters with a pass frequency band of 2 kHz to 30 kHz. This gave a total roll-off of 48 
dB per octave, w r hich was necessary due to large amounts of interference at 1800 Hz 
(caused by ventilation fans) and 31.25 kHz (produced by an unknown source in the 
building). 

7. Amplification and Filters - Later Work 

The later work was conducted during a portion of the year that was less noisy 
and included an isolation transformer for the filters which minimized the 60 Hz venti- 
lation fan noise. The PAR 113 Pre-Amplifier was used at a gain of 2000 for all drop 
diameters. The signal was then fed to a single Krohn-Hite 3202 R band pass filter, 
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passing frequencies between 1 Hz - 300 kHz for impact signal acquisition and 1 kHz to 
30 kHz for bubble signal acquisition. These band pass frequency ranges contain at least 
99 % of the total signal energy as determined by conducting trials both with and without 
filtering. 

8. Data Acquisition 

A digital data acquisition card, Computerscope, marketed by RC Electronics, 
was mounted in an IBM PC/XT and was used for data acquisition in all acoustical ex- 
periments. It is capable of sampling frequencies of up to 1 MHz at an amplitude resol- 
ution of 12 bits. It could sample up to 16 channels simultaneously, but the maximum 
sampling frequency would be lowered. The temporal resolution given by the 1 MHz 
sampling frequency was helpful when small time differences were measured from isolated 
impacts. A lower frequency of 250 kHz was selected for all of the bubble data samples, 
since a longer record length was needed than could be achieved at higher sample fre- 
quencies. 

9. Distance Measurements 

The early work used a video camera / tv monitor setup to measure the hori- 
zontal range from the impact to the hydrophone with an accuracy of + 1.5 cm. In the 
later work, distance was measured directly with a ruler, leading to an accuracy of ±0.5 
cm. 

10. Temperature Measurements 

Several Navy-issue thermometers were used in tins work with an overall range 
of -20 C to +50 C. Accuracy was ± 0.5 C. 

1 1. Salinity Measurements 

The salinity content of the water used for the saline dependence experiment was 
measured by an AGE Model 2100 Salinometer with an accuracy of 0.05 ppt. 

B. PHOTOGRAPHIC EXPERIMENT 

1. High Speed Photography 

Motion pictures of terminal speed 4.6 mm drops impacting the water surface 
were taken in conjunction with both this work and Snyder's (1990) work. A 400 frame 
per second Milliken camera was used to capture the time sequence from the initial im- 
pact through the violent tearing of the canopy. A complete analysis of this experiment 
is given in Snyder (1990) and is reviewed in this work in Chapter IV. 
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2. Video Camera Photography 

In the crossover energy experiment, it was necessary to see the cavity formation 
sequence for several velocities of a given drop size. This was accomplished with a Sony 
CCD-V99 video camera recorder. A 650 W movie light was used to provide adequate 
light to allow the use of a 1/1000 s shutter speed. This high speed was necessary to stop 
the motion of the crater formation. A 20 cm depth of field was possible when using the 
'macro' mode of the lens. The setup used is shown in Figure 6. 




Figure 6. Setup Used for the Crossover Energy Determination 
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III. NUMERICAL ANALYSIS METHODS 



A. INITIAL CALCULATIONS 

1. Signal Input 

The recorded acoustic energy was sampled at the predetermined time interval 
using the Computerscope data acquisition card. The digitized input has an amplitude 
resolution of 12 bits and is recorded as a voltage level which ranges between ±10 V. 
This voltage level is proportional to the pressure at a range determined by the geometry 
of the hydrophone and drop impact position. 

2. Correction to 1 m on Axis 

Both the impact and bubble signals were then corrected to an equivalent pres- 
sure which would be at 1.0 m and on the vertical axis below the drop impact point. This 
correction must take into account the 1/r spherical spreading as well as a cos 6 depend- 
ence as shown in Figure 7. This latter correction assumes that both the bubble and 
impact behave as acoustic dipoles. The bubble radiation pattern has been confirmed for 
smaller, terminal velocity drops at normal incidence (Kurgan, 1989). If the original 
voltage amplitude is given by A [v], then the correction factor is: 



A 



at\m on axis 



A x 



h 2 + z 2 
100 z 



( 4 ) 



where z is the depth of the hydrophone in cm and h is the horizontal distance from the 
hydrophone to the drop impact point in cm. 
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Figure 7. Geometry for Dipole Correction for Range and Angle 
3. Near Field Correction 

The near Field effect caused by being within several wavelengths of a dipole 
causes an error between the measured pressure field and the extrapolated (at 1 m) pres- 
sure field. The proper correction factor has been derived to be (Medw’in and Beaky, 
1989): 



^at lm on axis far field ^ at 1 m on axis 



where k is the wave number in m~ l and R is the range from hydrophone to impact in 
m. In further equations, a prime superscript will be used to denote "at 1 m on Axis, Far 
Field." 

The near field correction was not made for the earlier work due to the 
hydrophone placement at a depth of 15 cm, which is at least one wavelength from the 
surface for all frequencies above 10 kHz. These measurements were taken in the 
redwood tank before the anechoic lining w r as installed, and the reverberation clearly 
dominated any errors introduced by the near field values. 

The later work required the near field correction since all frequencies of interest 
were in the near field at a hydrophone depth of 6 cm. 
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4. Fourier Transform 

The temporal signal was converted to its frequency components by a 2048 point 
Fast Fourier Transform. The routine used (from the Borland Turbo Pascal Numerical 
Toolbox) was based on the Coolev-Tukev method and is given by: 



X(k) = A 7- 1/2 e ~ a ™ kls 

A= 0 



( 6 ) 



Some methods of FFT do not include the J~N in the forward transform and use a \ jN 

term in the inverse transform instead. Also, a factor of may differ between this 

2n 

method and others. These factors are the most likely reason why results obtained from 
another source may differ slightly with the results presented in this paper. 

The unit analysis continues: 



A' 
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( 7 ) 



where bin width is given by 



Vm\Vidih[Hi\ = 
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The sensitivity of the hydrophone is uniform over the frequency range of interest and 
has a value of -91.5 dB re IV/ Pa. This is equivalent to 37,580 Pa'V. 

With this information, the desired result of spectral density for a 1 Hz band- 
width at 1 m on axis in the far field can be obtained. For simplicity, these conditions 
will be denoted by a double prime superscript. For a case with a gain of 2000, 
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The result of this calculation is the spectral density from a single bubble or impact. 

B. DATA ANALYSIS PROGRAMS 
1. Single Drop Processing 

All of the initial calculations were carried out by a single program written to run 
in batch mode (no user interaction). Each bubble and impact were processed sequen- 
tially and yielded an output of frequency versus spectral level [Pa iJJFz ]. Each individual 
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spectrum was evaluated for noise levels, peak frequency, and the presence of secondary 
bubbles and atypical signals. The duration of the bubble record was fixed at 8.2 ms to 
establish a standard record length for subsequent processing and was based on the 
shortest record length able to capture at least 99% of the energy from the bubbles 
produced by the 4.6 mm drops. The impact record length was set at 256 microseconds, 
long enough to capture the impact for all drop diameters. 

2. Multiple Drop Processing 

The spectral densities of a number of bubbles or impacts with a common pa- 
rameter (drop diameter, temperature or salinity) were averaged by a second, user-written 
program. The output of this program is frequency versus spectral density [Pa 2 / Hz], 
which can then be plotted or used for energy calculations. The reasons for requiring an 
average value are discussed below. 

3. Plotting 

The results from the averaging program were plotted using Quattro Pro, which 
was selected for its versatility in vector (column) calculations and rich graphic capabili- 
ties. For each set of data, the total spectral density is calculated by adding the spectral 
density from the impact together with a fraction of the spectral density from the bubble. 
This ratio takes into account that impact energy is produced with each raindrop strike, 
while bubbles are created and radiate energy for only a fraction of the drops. 

C. DATA AVERAGING 

Due to the natural variability in sound spectrum shapes produced by impacts and 
bubbles, it is necessary to view the results using statistical techniques. No raindrop will 
produce a spectral shape exactly like another. It is therefore necessary to average the 
bubble and impact spectrums in both frequency and ensemble. The frequency average 
consists of a 1 kHz wide moving filter applied to smooth the individual spectrum before 
plotting. The ensemble average uses all drops of a similar category (e.g., all 4.6 mm 
drops). Also, the standard deviation of the spectral levels at the dominant bubble fre- 
quency is given, as well as the standard deviation of the total energy per raindrop, to 
present the reader with a good measure of the variability of energy present in rainfall. 

D. ADDITIONAL CALCULATIONS 

1. Average Spectral Energy per Raindrop 

To obtain the spectral energy per raindrop, the frequency dependant average 
spectral density [/V/Hz] at 1 m on axis due to a given diameter of raindrop (or other 
parameter) is integrated with respect to frequency. The integration will yield the square 
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of the pressure radiated by that average raindrop. Assuming a dipole radiation pattern, 
this can then be converted to spectral intensity on axis and integrated by using ring ele- 
ments [(2 ;rR sin 6) R d0] over the area of radiation. It follows that 



EnergyiJ) 
Drop Time 




A" cos 2 d sin ddd 



which, for R= 1 m yields 
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where Time is given by the length of time the bubble radiates sound. For this analysis, 
the time will be the record length of data extracted to meet the requirements of the FFT 
algorithm (# of data points is a power of two) and has been chosen to be 8.2 ms. A 
record length of 8.2 ms has been shown in the laboratory to capture at least 99% of the 
total energy radiated by a typical bubble. 

2. Calculated Underwater Sound Spectrum Levels due to Rainfall 

For a known drop size distribution, it is possible to use the average spectral 
density at lm on axis per raindrop computed previously to obtain the calculated rainfall 
spectrum, RS(f), for the given drop size distribution. This is a straightforward calcu- 
lation: 



RS(f) = ^ a " rr ( d ) T s 

D 



(Effective R) 2 



( 12 ) 



where RS(f) has units of Pa 1 ! Hz, A" has units of Pa 1 ! Hz jdrop, RR(D) is the rainfall rate 
dropsjm 1 s of diameter D and T, is the total sample time used for data collection (8.2 
ms for this experiment). The effective R is determined by the location of the raindrop 
with respect to the hydrophone, taking into account the dipole radiation pattern of the 
raindrop. This calculation allows one to compare the theoretical results with previously 
measured spectrum levels if the hydrophone geometry is known. An example will be 
shown in Chapter VI. 
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IV. DESCRIPTION OF RAINFALL MECHANISMS WHICH CAUSE 

BUBBLE FORMATION 



A. TYPE I MECHANISM 
1. Background 

Bubbles generated by the Type I raindrop mechanism have been previously de- 
fined (Snyder, 1990) as those bubbles caused by terminal velocity drops sized from 0.8 
mm to 1.1 mm diameter. Bubbles of this mechanism are the primary component of low 
rainfall rate and occur during light wind and at near normal incidence (Nystuen and 
Farmer, 1987). They produce a spectrum with a sharp peak centered at about 15 kHz. 
This peak has been shown to broaden and weaken as wind speed increases, due to a 
lower production rate of bubbles when drops impact at oblique incidence (Kurgan, 1989; 
Medwin, et al.,1990; Nystuen, 1991). 

The range of terminal velocity drops that produce Type I bubbles is shown in 
Figure 8. This figure shows all of the previously identified regions of bubble production. 
The curves labeled "NCPA" have been described by the National Center for Physical 
Acoustics in Mississippi (Pumphrey and Elmore, 1990). Most of the work done by 
NCPA has been away from the terminal velocity curve. Since rain always falls at ter- 
minal velocity, research conducted away from the terminal velocity curve does not rep- 
resent rainfall at sea. 
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Figure 8. Known Bubble Production Regions 
2. Bubble Production Mechanism 

Type I bubbles caused by drops falling at terminal velocity have been thor- 
oughly researched both experimentally (Kurgan, 1989; Pumphrey, et. al., 1989; Meduin, 
et.al. 1990) and theoretically (Oguz and Prosperetti, 1990; Longuet-Higgens, 1990). The 
exact mechanism of bubble production is not known, but several theories present likcl\ 
conclusions. 

After the impact of the drop, a conical cavity is formed. The shape of the cavity 
is the key visual identification of the Type 1 mechanism. For a small range of drop sizes 
at terminal velocity (0.8 - 1.1 mm), a bubble is pinched off at the bottom of the cone 
after a variable delay time of about 20 ms. 1 his pinch-ofT may occur because of capillan 
waves that join at the bottom of the conical surface (Longuet-Higgens, 1989). This 
theory states that the pinch-off will occur when the conical angle reaches 109.5 degrees. 
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Another possible explanation is that the bubble is pinched ofT when the conical cavity 
forms a "nipple" at the bottom (Oguz and Prosperetti, 1989). Under certain geometrical 
conditions, this "nipple" separates from the cavity to form a small bubble. 

Both of these theories require a delicate geometrical structure to allow creation 
of a bubble. However, when a drop strikes a flat surface at oblique incidence, which is 
equivalent to a drop falling against a wave at other than normal incidence, the percent- 
age of bubbles produced decreases rapidly (Kurgan, 1989; Mcdwin, et. al., 1990). All 
theories support the conclusion that a bubble is pinched off only when symmetrical 
conditions apply. Results from rainfall studies also show that when wind speed increases 
(causing oblique incident angles), the characteristic 15 kHz peak created by Type 1 
bubbles decreases in amplitude and becomes broader in the frequency domain. 

The Type 1 mechanism (pinch-off) has also been observed in drops at velocities 
other than terminal velocity (Pumphrey, 1989). These drops are described by the "reg- 
ular entrainment" region of Figure 8 and exhibit a similar 15 kHz peak. Drops in this 
region do produce bubbles 100% of the time but are not useful in the study of rainfall, 
since rain is comprised of only terminal velocity drops. However, bubbles produced by 
drops in this region may be commonly present in the ocean due to the spray from waves 
or as secondary bubbles during heavier rainfall. 

B. TYPE II MECHANISM 
1. Background 

Bubbles caused by the Type II mechanism have been previously defined 
(Snyder, 1990) as those dominant bubbles caused by terminal velocity drops ranging in 
size from 2.4 mm to 4.6 mm diameter (Figure 8). This paper broadens the definition to 
include terminal velocity drops between 2.2 mm and 4.6 mm and secondary bubbles as 
well. Bubbles produced by this mechanism were first studied by Franz who stated, "The 
bubble component of the underwater sound energy from the splash of a water droplet 
is very erratic under most conditions." (Franz, 1959) This notion of unpredictability was 
carried on by Pumphrey and Crum in 1989 when this mechanism was labeled "irregular 
entrainment." 

Recent work has shown that although a wide variety of drop size / impact ve- 
locity combinations will form bubbles by the Type II mechanism, statistically predictable 
results can be obtained when terminal velocity drops (rainfall) are investigated (Snyder, 
1990). 
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Results presented later show that formation of a dominant bubble varies from 0% at 
the low end of drop sizes researched (2.2 mm diameter) to a peak of 65% at approxi- 
mately 4.2 mm diameter. Drops above 4.6 mm diameter were not considered, as Nature 
rarely produces drops larger than these. The frequency of the dominant bubble 
produced by this mechanism can be directly related to drop size (Snyder, 1990), which 
indicates that 2.2 mm to 4.6 mm drops falling at terminal velocity are good candidates 
for estimating the distribution of raindrop sizes by acoustic means. A typical time do- 
main sequence of a dominant bubble created by a Type II mechanism is shown in 
Figure 9. 




Figure 9. Time Domain Portrait of a Dominant Bubble 
2. Bubble Production Mechanism 

Although Franz had photographed the Type II mechanism in 1959, only recent 
photographic efforts described by Snyder have identified the sound producing compo- 
nents of this mechanism. Sketches from the high speed photography work are shown in 



20 




Figure 10 for a 4.7 mm drop (Snyder, 1990). The time between frames is 2.5 ms, al- 
though the frames shown are not necessarily consecutive. 

The first frame shows the drop slightly prior to impact. Note that the drop is 
not spherical but is flattened on the bottom. This shape has been predicted theoretically 
(Pruppacher and Fitter, 1971). 

The next frame shows the formation of the crown and the beginning of a 
hemispherical crater in the water, which begins shortly (2 ms) after impact. Note the 
spray of aerosols ejected by the upward travelling water mass. The next several frames 
show the closing of the canopy, most likely due to surface tension. The canopy closes 
about 15 ms after the impact. 

After the canopy closes, a jet of water continues to rise from the site at which 
the canopy closed. The next frames show the growth of the canopy to a height of 2.9 
± 0.3 cm at a time of 25 - 30 ms after impact. When the jet has reached a peak in height, 
a downward moving jet appears at the base of the upward moving jet. This jet plunges 
downward towards the bottom of the flattened crater with an observable cant (angle). 

The next frame shows the jet piercing the bottom of the fiat-bottomed crater. 
Several ms later, a bubble is broken off from the tip of the turbulent jet, apparently due 
to buoyant forces that act on the air trapped in the tip of the jet. Delay times from the 
impact to the onset of the bubble range from 35 ms to 65 ms and depend on drop di- 
ameter (Figure 11). 

The jet then retracts into the crater, which is already in the process of collapse. 
The canopy remains on the surface for 20 - 50 ms before tearing itself apart, violently 
spraying aerosols as far as 5 cm from the point of impact. 

From the photographic evidence seen, the necessary criterion for bubble for- 
mation is the cant of the turbulent jet as it plummets downward through the crater. If 
the jet is perpendicular to the water surface, no buoyant forces are exerted, and no 
bubble is produced. It is believed that drops incident at oblique angles would have a 
higher percentage of bubble production. From a simple geometric standpoint, it can be 
argued that oblique incidence would provide the imbalance that is necessary for bubble 
production. Wind, therefore, should enhance the spectrum produced by this Type II 
drop mechanism. This is in direct contrast to the Type I mechanism, which shows a 
decrease in bubble production in the presence of wind. 
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Figure 10. High Speed Photographic Sequence of Type II Mechanism 
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Figure 11. Delay Time From Impact to Onset of Bubble 



3. Secondary Bubble Production 

In a small number of cases, a second, smaller bubble is produced by the Type 
II mechanism. A time domain example of this is shown in Figure 12. This phenomenon 
deserves attention because of the energy potentially contributed by the secondary' bub- 
bles to the overall spectrum. 
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Figure 12. Time Domain Portrait of a Secondary Bubble Which Precedes a Domi- 
nant Bubble 

4. Wobble Production 

In a moderate number of cases (more fully described in the next chapter) a 
strange phenomenon resembling a gated sinusoid appears (Figure 13). This phenome- 
non has been labeled a "wobble" because of its appearance. The mechanism of sound 
production has not been identified, but it may be related to pulsations of the turbulent 
jet or caused by capillary waves modulating the jet. The wobble has been observed both 
before and during the dominant (as opposed to secondary) bubble. On several occasions 
the wobble has been the only source of sound. 

The start of the wobble is very abrupt, as is its end. The start has been observed 
as early as 16 ms after impact, and as late as 1 ms prior to the dominant bubble (the 
dominant bubble delay time is a function of drop diameter). It maintains an almost 
constant amplitude during the entire event, which can last between 5-13 ms, and has 
approximately the same duration as the dominant bubble. The amplitude, however, re- 
mains constant and does not show an exponential decay, as does the bubble. The 
wobble has been investigated to determine the amount of energy contributed (Chapter 
V). In addition, the frequency of the wobble and its entire spectrum has been examined 
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to see whether one can separate the wobble from the dominant bubble in the frequency 
domain (Chapter V). 




Figure 13. Time Domain Portrait of a Wobble (starts at 43.4 ms and continues into 
the dominant bubble) 

C. CROSSOVER KINETIC ENERGY 

1. Background 

Franz theorized that the difference between the two mechanisms that he had 
seen (non-terminal velocity experiments which we assume were Type I and Type II) was 
based on drop velocity. It is easy to see from other work, such as the photographs in 
Pumphrey's (1989) dissertation, that velocity is not the only criterion. This is also evi- 
dent from the regions identified by NCPA (Figure 8). A possible alternative to velocity 
is the drop's kinetic energy at impact. 

2. Experimental Setup 

To test the theory that kinetic energy is the critical parameter for type II bubble 
creation, several readily available points of data were used. Our research has shown that 
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a terminal velocity 2.2 mm drop produced a bubble from the Type II mechanism 0% 
of the time. Drops larger than 2.2 mm diameter falling at terminal velocity often 
produced bubbles. It was therefore assumed that the 2.2 mm drop possessed the mini- 
mum kinetic energy necessary to initiate the Type II mechanism. 

Our postulate was that if a drop larger than 2.2 mm, falling at a lower velocity, 
had the same (or less) kinetic energy as the 2.2 mm drop at terminal velocity, the Type 
II mechanism would not be present. If the drop were released from a higher point 
(giving a larger kinetic energy) the Type II mechanism would be seen. If the drop were 
released from a lower point, the Type I mechanism would be seen (see Figure 8). 

The 5.2 mm drop diameter was selected for this experiment because of the 
availability of eyedroppers that yielded this drop size. 100 drops were released from 2.15 
m and had a calculated velocity of 5.4 m's at impact (corresponding to the XCPA "ir- 
regular entrainment mechanism"), using the formula (Pumphrey, 1989): 



where v T is the terminal velocity of the drop. Drops were also released from 0.25 m 
(corresponding to the lou'er kinetic energy of the Type I "regular entrainment region") 
to check for the Type I mechanism. These lower velocity drops had a calculated impact 
velocity of 2.2 m/s. A video camera was used to provide a record of the experiment as 
well as to provide frame-by-frame playback capability to study the drops. 

3. Results 

At 2.15 m height of release, several of the characteristic Type II mechanisms 
were noted for the 5.2 mm drop (velocity 5.4 m/s). The hemispherical crater and crown 
were easily seen. The remainder of the drops formed another structure, unlike Type I 
or Type II, which did not appear to produce bubbles. This type of drop splash has been 
noted previously by Franz and w r as identified as the "third type" (Franz, 1959). 

At a height of 0.25 m the results for the 5.2 mm drop (velocity 2.2 m/s) were 
quite different. Most of the drops had the characteristic conical shape of the Type I 
mechanism, and no Type II mechanisms were noted. Several of the drops had the "third 
type," which may be an intermediate structure between Type I and Type II. 

However, further investigation yielded similar results for the 5.2 mm drop at a 
height of 1.25 m (corresponding to an impact velocity of 4.6 m/s). From Figure 8, only 
the Type II mechanism should be present at this velocity. The corresponding kinetic 
energy is well above that of the 2.2 mm drop at terminal velocity. 




(13) 
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These results appear to support the claim made by the NCPA that the region 
labeled "irregular entrainment" in Figure 8 is indeed properly labeled for non-terminal 
velocity drops. There does not seem to be a clear, unique dividing line of kinetic energy 
separating the Type I and Type II mechanisms for all drop diameters and all velocities. 
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V. IMPACT AND BUBBLE ENERGY ANALYSIS 



A. ENERGY SOURCES 

The energy radiated into the water by a drop splash consists of four primary ele- 
ments: Impact sound, underwater hydrodynamic motion, surface wave generation and 
bubble sound radiation. Only two components, impact and bubble sound radiation, are 
of present interest. Hydrodynamic effects, however, can contribute substantially to the 
overall near field pressure, especially for the impact water hammer followed by the gen- 
eration of the crater. The subject of impact hydrodynamics continues to be researched 
at this laboratory. 

B. SAMPLE SPACE 

A total of 202 drops which produced bubbles and the corresponding impacts were 
analyzed. Drop temperature and surface temperature were at 20 °C (room temperature) 
and fresh water was used. The distribution of drop sizes is given in Table 1. The 2.6 
mm drop data was not used in further analysis due to the small sample space. Other 
drop sizes (except the 3.7 mm diameter) had a standard deviation of peak bubble fre- 
quency (frequency at which the spectral density peaked) of less than 500 Hz (Snyder, 
1990). The 3.7 mm diameter had a standard deviation of 1 kHz. The selection of ap- 
proximately 30 drops of each size was based on the desire to obtain a large enough 
sample space to estimate the standard deviation of the dominant bubble frequency. The 
selection criterion was based on data obtained for previous work (Snyder, 1990) and is 
adequate for this work. 

Drop sizes were fixed by the availability of eye droppers and the desire to obtain a 
difference of about 0.3 mm between each drop diameter. 

C. IMPACT ENERGY 

1. Impacts of Small (0.8 - 1.1 mm) Drops 

The impact would be expected to contribute to the overall energy radiated by a 
drop, since the impact noise is present for each drop strike. The impact energy of small 
(1 nun) raindrops has been previously investigated (Kurgan, 1989). As expected, the 
radiated impact energy from this size of drop was small, on the order of 0.01 picoJoules. 
Kurgan's impact spectrum for small drops peaked at approximately 15 kHz and had a 
-3 dB bandwidth of about 14 kHz (Figure 14). The wide spectrum was expected, as the 
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Table 1. DROP DIAMETER SAMPLE SPACE 



Drop Diameter (mm) 


Number of Drops 


2.6 


3 


2.7 


29 


3.1 


36 


3.4 


33 


3.7 


14 


4.0 


28 


4.2 


29 


4.6 


30 



short duration of the impact (typically 50 - 100 /rsec) in the time domain translates into 
a broad spectrum in the frequency domain. However, the appearance of the peak in 
Kurgan's work is due to the impact spectrum not being corrected for the high and low 
pass filters used (as noted by Kurgan). When corrected, the spectrum has a trough at 
about 6 kHz and then has an increasing magnitude with with decreasing frequency, as 
opposed to that shown in Figure 14. This characteristic of the impact has been de- 
scribed by Pumphrey (1991) and has been experimentally verified for other drop diam- 
eters in this laboratory. 
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